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Abstract 

We measure the mass of the rj meson using iJj{2S)^r]J/^j events acquired with the CLEO-c 
detector operating at the CESR e+e" coUider. Using the four decay modes >77, Sir^, n~^Tr~Tr^, 
and 7r+7r~7, we find M^=547.785±0.017±0.057 MeV, in which the first uncertainty is statistical 
and the second systematic. This result has an uncertainty comparable to the two most precise 
previous measurements and is consistent with that of NA48, but is inconsistent at the level of 6.5a 
with the much smaller mass obtained by GEM. 
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The T] meson, the second-hghtest pseudoscalar, is commonly understood as being predom- 
inantly in the SU(3)-fiavor octet with a small singlet admixture, so that it has comparable 
uu, dd, and ss content and virtually no gluonium component [H, [1]. Its mass is of fundamen- 
tal importance to understanding the octet-singlet mixing as well as the gluonium content 
of both 77 and rj' [sf, although theoretical and phenomenological precision on related pre- 
dictions [1, has not yet matched that of experiment. On the experimental side, there 
has long been a situation of conflicting measurements that improvements in precision 
have been unable to resolve. Indeed, the current status is the worst it has ever been, with a 
confidence level (CL.) of 0.1% that the measurements are consistent j^. In 2002 the world 
average was M^=547.30±0.12 MeV Q and results included in it were generally consistent 
with one another, but only because 1960's-era bubble chamber experiments, all of which 
favored a larger M^, had been dropped Then NA48, based on exclusively reconstructed 
?7^37r° decays, reported M^=547.843±0.030±0.041 MeV Q, which appeared to vindicate 
the dropped experiments. In 2005 GEM reported M^=547.311±0.028±0.032 MeV [ToJ us- 
ing the ^He recoil mass in p+d^^B.e+X . The GEM result was consistent with less precise 
results made during the period 1974-1995, but also was eight standard deviations below that 
of NA48. More measurements with sub-100 keV precision are needed to clarify the matter. 

This work presents a new measurement of using if) [28)^7] J /ip. Events were acquired 
with the CLEO-c detector at the CESR (Cornell Electron Storage Ring) symmetric e^e~ 
collider. The data sample corresponds to ~27 million produced ip{2S) mesons, of which 
about 0.8x10^ decay to rjj/ip. We measure the mass by exploiting kinematic constraints 
in the decay chain ip{2S)—^rjJ/ip, J/ip^C^l' (£='==e^ or /x^), and 77— i>77, 37r°, tt+tt^tt^, or 
7r+7r~7. Because both ip{2S) and J/ip are very narrow resonances with precisely known 
masses, the constraints enable a significant improvement in t] mass resolution over that 
achieved by the detector alone. This is the first measurement to use tlj{2S)^rjJ/ijj. 



The CLEO-c detector is described in detail elsewhere [ll[; it offers 93% solid angle 
coverage of precision charged particle tracking and an electromagnetic calorimeter com- 
prised of CsI(Tl) crystals. The tracking system enables momentum measurements for par- 
ticles with momentum transverse to the beam exceeding 50 MeV/c, and achieves resolution 
ap/p ~0.6% at p=l GeV/c. The barrel calorimeter reliably measures photon showers down 
to ^^=30 MeV and has a resolution of a^/E ^5% at 100 MeV and 2.2% at 1 GeV. 

Event selection begins with that described in Ref. [13] for the four rj decay modes used 
here. We also accumulate samples of ip{2S)^n^7i~J/ilj, 7r°7r°J/?/', and tt'^J/i/j for studies 
of systematic uncertainties. Every particle in the decay chain is sought, and events are 
separated into those with J /ip—^ji'^ and J /ip^e'^e~ . Leptons are loosely identified and 
restricted to | cos^^±|<0.83, where 9 is the angle of the track with respect to the incoming 
positron beam. Lepton momenta are augmented with calorimeter showers found within a 
100 mrad cone of the initial track direction, under the assumption that they are produced 
by bremsstrahlung. All photon candidates are required to be located in the central portion 
of the barrel calorimeter where the amount of material traversed is smallest and therefore 
energy resolution is best: | cos6'^|<0.75. Backgrounds of 1-4% are present in each rjj /ip sub- 
sample [l2|, consisting of cross-feed from other rj decays as well as other ip{2S) — * XJ/%Ij 
decays: tt+tt'J/t/;, 7r°7r°J/?/', and 7Xcj, Xcj^lJ/'^l^- 

Kinematic constraints are applied in two two-step fits: first, the lepton tracks are con- 
strained to a common origination point (vertex) and thence to the J /if) mass, MjJ^ j6[; 
second, the constrained J/ip, the beam spot and the 7? decay products are constrained to a 
common vertex and then to the il>{2S) mass, M^^^s) pl- Separate fit quality restrictions are 
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applied to vertex {xD mass (Xm) constraints. For rj^n'^n'Ti'^, the 7r°— i>77 candidate is 
constrained to the 7r° mass prior to the fits described above. The decay tj^Stt^ is treated as 
77— i>67 because rehably making a unique set of correct photon- 7r° assignments is not possible; 
typically several such assignments per event of comparable probability exist and are indis- 
tinguishable. To ensure that only the best measured events survive into the final sample, 
the restrictions from Ref. [l2| are tightened to Xv/d-o.f.<10 and Xm/d-o.f.<5 for both 
the J/ip and ip{2S) constrained fits. 

Alternative event topologies are used to compare measurements of the 7r° mass to 
its established value, M^o°*^=134.9766±0.0006 MeV for two different 7r° momentum 
ranges. The first (vr^) is tlj{2S)-^TC^J/tlj, 7r°— i>77, which features a monochromatic vr^ with 
p~500 MeV/c, and the second (vr^) ilj{2S)^rjJ /ip, 77— >-7r"'"7r~7r*', 7r°— ^77, which contains vr^'s 
with p~ 0-250 MeV/c. For these tests, the individual photons (instead of a constrained vr*^) 
are used in the M^{2S) constraint on all final state four-momenta. 

Each event yields an invariant mass M of the kinematically-constrained decay products; 
a single mass value is extracted for each decay mode i by fitting a Gaussian shape to the 
distribution of ^j=Mj— Mq, where Mq is simply a reference value, either the current Particle 
Data Group world-average M^^*^=547.51±0.18 MeV j^, or, in the case of the 7r° cross- 
check modes, M™^ as given above. The fits are restricted to the central portion of each 
5 distribution because the tails outside this region are not represented well by a single 
Gaussian form. The fits span ±1.6o" to ±2. Oct about the peak (5), where a is the fitted 
Gaussian width, and in all cases the resulting fit has a C.L. exceeding 1%. The distributions 
of 6i for 7r° and 77 decay with fits are shown in Figs. [1] and [21 respectively. Other shapes that 
might fit the tails, such as a double Gaussian, have been found to yield unstable fits and/or 
do not improve precision of finding the peak. 

There is an unavoidable low-side tail in any monochromatic photon energy distribution 
from the CLEO calorimeter. It originates from losses sustained in interactions prior to 
impinging upon the calorimeter and from leakage outside those crystals used in the shower 
reconstruction. This asymmetric photon energy resolution function also results in a small 
but significant systematic bias in (5): for simplicity of the kinematic fitting formalism, input 
uncertainties are assumed to be symmetric, and a bias occurs if they are not. This bias in 
fitted Gaussian mean is mode- dependent because each presents a different mix of charged 
and neutral particles. 

The biases I3i are estimated by following the above-described procedure on MC signal 
samples. Each (3i is the difference between the Gaussian peak value of the distribution 
and the input M^^. We define the bias as [3i={5i)y[c^ in which we use the MC input 
instead of M^^^ for Mq. A non-zero value of [3i means that, for mode z, the Gaussian peak 
mass is offset from the true mass and must be corrected. We evaluate the four /3j for 
trial values of M^'^ (547.0, 547.3, 547.8, and 548.2 MeV) that cover the spread of previous 
measurements. The biases extracted for these inputs are consistent, and the final bias 

for each mode, shown in Table [H is taken as their average. Bias values for the 7r° and 7r° 
cross-checks are determined similarly. 

Table [T] summarizes results by decay mode. Both 7r° and 7r|^ mass values are consistent 
with expectations within their respective statistical uncertainties. The total number of 
reconstructed events involved in the determination of is 16325. The four values of 
— Pi have an average, weighted by statistical errors only, of {{S) — /?)s=277±17 keV with 
a x^=4.8 for three degrees of freedom (C.L.~20%). 

Systematic errors are summarized in Table [ITl Uncertainties that are uncorrelated mode- 
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to-mode, including statistical, are used to determine the weights (wj=0.19, 0.03, 0.60, and 
0.18 for 77, Stt", 7r"'"7r~7r°, and 7r"'"7r~7, respectively) applied to combine values from the four 
modes into the weighted sum ((5)-/?)w=Ef=i^^iX((^i)-A)=272±17 keV. 

As the mass distributions are not perfectly Gaussian, there is some systematic variation 
of the peak value with the choice of mass limits for each fit. However, as long as the limits 
chosen do not result in a confidence level below 1% and remain roughly symmetric about the 
peak, such variation is observed to be bounded by approximately half of a statistical stan- 
dard deviation. Hence this value was assigned as a conservative estimate of the systematic 
uncertainty attributable to the fit limits. 

Uncertainties attributable to imprecision in the masses of the J /if) (11 keV) and ip{2S) 
(34 keV) mesons ^ are directly calculated by repeating the analysis using an ahered ip{2S) 
or J/ip mass and the deviation in (6) per "la" change from nominal taken as the error. 

The bias /3j from kinematic fitting can be attributed to the effect of the asymmetric 
resolution function of photons; the bias for the four modes varies by an order of magnitude, 
with larger values corresponding to modes with more photons. The 7r° cross-check indicates, 
within its 47 keV statistical precision on {6)— (3, that the bias is indeed accurately estimated. 
A more sensitive cross-check comes from comparing the shifts in data and MC for 
7] 71^71^71^ when the vr'' mass constraint is removed: the MC predicts an increase in bias 
of 76±6 keV compared to an observed shift in the data of 55±24 keV {i.e. the shift in 
bias of 76 keV is verified in the data within the 24 keV statistical error, which amounts to 
about a third of MC shift itself). Based on these comparisons and the generally favorable 
agreement [13] between data and MC characteristics, we take one third of the bias central 
value {Pi/ 3) as our estimate of the systematic uncertainty. 

By performing mass fits on MC signal samples with and without simulated backgrounds, 
it is determined that for mode i, modeled backgrounds reduce the bias by the amount Bf. 
21±19, -3±67, 2±17, and -13±27 keV for the 77, 3n°, n+n-rr^, and tt+ti-j channels, 
respectively, where uncertainties listed are statistical. The unmodeled background in the 
7r+7r~7 sample jl2| is estimated to have an effect on that is negligible compared to the 
27 keV uncertainty on the modeled background for this mode. After weights are applied, 
the net effect is a positive offset to M^, of B^=3±12 keV. 

Uncertainties in charged particle momentum and calorimeter energy scale are evalu- 
ated by shifting those scales by the appropriate amount and repeating the analysis. The 
charged particle momentum scale is confirmed at high momentum (~1.5 GeV/c) using 
unconstrained J/ip^n'^fi~ decays from ^{2S)^7i~^Tr~ J/ip and J/ip events. A low- 
momentum (75-500 MeV/c) calibration, which is more relevant to our measurement of 
M^, can be made by comparing the mass of ip{2S)~^7i^7i^ J/ip to M™^-) with no kine- 
matic constraints on the Tr"*"??" but with a mass-constrained J/il!^i~^i~; this checks the 
TT+TT" momentum scale because the accurately known J/ip mass takes up ~84% of the 
available energy. Events of both types are subjected to Gaussian fits to mass difference 
variables similar to 6, in the first case to 6^+fj,-=M{^'^ ^~)—MjJ^ , and in the second to 
(5^+^-=M(7r+7r~ J/-?/;)— M^^'^^-j. For MC simulation, where we can employ our perfect knowl- 
edge of the magnetic field for the momentum scale, the means and statistical errors of these 
measures are (5^+^-)=— 90±22 keV and {6.„+-„-)=2±3 keV. This demonstrates that this 
technique is accurate to |(5^+^-)|/Mj/^~3xlO~^ and | (^tt+tt-) |/(^v(25) ~ ^j/V')— 1^ 10~^- 
The magnetic field scale in data is tuned to that value which keeps both of these means close 
to zero: with this setting, measurements yield (5^+^-)=7±46 keV and {6t,+t,-)=—23±6 keV, 
indicating a similar level of sensitivity as the MC samples. Therefore we quote a relative 
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momentum scale accuracy of 0.01% and use this value for our la systematic variation. 

Several processes are used for the calorimeter calibration: inclusive decays to 77 (where 
we can constrain M(77) to a known mass), e'^e~ (in which energy-momentum 

conservation and well-measured track momenta allow constraint of the photon energy), and 
ilj{2S)^'~fXcj (where the transition photon energies are known well). The photon energy scale 
is calibrated to give the correct peak, i.e. the most probable value, for any monochromatic 
photon energy distribution. These calibrations are combined and result in an overall energy 
scale known to 0.6% or better over the energy range (30-400 MeV) relevant for photons from 
the slow r^-mesons produced in ilj{2S)—>-r]J/tlj. 

Any deviation from ideal in momentum or energy scale is substantially damped by the 
mass constraints, as is evident from Table [Tll the relative momentum (energy) scale uncer- 
tainty of 0.01% (0.6%) induces, at most, ~1 (~5) parts in 10^ shift in //-mass scale. 

We have also computed {{6) — /5)s when the decays occur in combination with either 
J/ilj—>-fi~^ fi~ or J/ilj—>-e^e~ separately; its value for J/tp^e~^e~ events is higher than that 
for J/ijj^ij,~^ fi~ by 98±34 keV (+2.9(t), where the error is statistical only. Broken down 
by mode, this difference is 71±57 keV (+1.2a) for 77, -82±208 keV (-0.4a) for Stt", 
119±49 keV (+2.4ct) for 7r+7r"7r°, and 141±91 keV (+1.5ct) for ir+n-'y. We do not observe 
such an effect in MC simulations. Further investigations, detailed below, reveal no firm 
explanation. To allow for a hidden systematic effect, we add an MC modeling uncertainty 
of 46 keV; it is the dominant uncertainty in this analysis. 

In order to investigate the possibility that there could be an unmodeled systematic 
pull of J /ip^e^e~ or J/-?/;— /i~ decays in the kinematic fitting process, we examine 
ip{2S)^-K^-K^ J /ip decays. If such an effect existed, we would expect to observe a shift 
between the dipion mass computed with the constrained pion (vr^) momenta relative to 
the unconstrained (vr^) values. However, we find no evidence for such a pull: defining 
/£+£- = (M(7r+7r~)— M(7r+7r~))f+£-, the difference fe+e- — //x+a(-=13±11 keV in data and 
7±9 keV in MC simulations (errors shown are statistical). 

To investigate the effect of less-well-measured events upon the analysis in general and 
the e~^e~-fi~^ fi~ discrepancy in particular, we have repeated the analysis after tightening 
the kinematic fitting restrictions from x^/d.o.f.<10 and Xm/^-^-^-<^ J/i^ ^(25') 
constrained fits to 5 and 2, respectively, losing about 40% of the original events. The net 
offset from backgrounds changes from 3±12 keV to — 8±14 keV. The overall final t] mass, 
including the background offset, changes by — 16±17 keV, demonstrating stability of the 
measured mass with respect to the kinematic fit quality. For this restrictive selection, the 
Mrj difference between e~^e~ and /i"*"//" events in data goes down to 48±44 keV, whereas the 
MC difference remains near zero. It appears that the e"'"e~-/x"'"/i~ discrepancy moderates for 
this class of events, but the statistical precision is not conclusive. 

In order to study dependence of upon the time of data collection, we divide the data 
into nine contiguous and consecutive data-taking periods. One mass from each period is 
obtained by averaging the results obtained from the four modes with statistical weights. 
The foi' the nine values to be consistent with their statistically-weighted average is 9.5 
for 8 degrees of freedom (C.L.=25%), demonstrating the absence of any time- dependent 
systematics. 

After combining the ((5i)— A) values in Table [T] using the quoted weights, including the 
aforementioned net effect of backgrounds B„, and adding the M^^*^ offset, our result is 
M^={{5)-(3)^+By,+M^^^= 547.785±0.017±0.057 MeV, where the first error is statistical 
and the second systematic. This result has comparable precision to both NA48 and GEM 
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measurements, but is consistent with the former and 6.5 standard deviations larger than the 
latter. All four prominent rj decay modes contribute to this result, and each independently 
verifies a significantly larger than obtained by GEM (~2.0(T in Stt", more for each of the 
other three) . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences 
and Engineering Research Council of Canada. 
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TABLE I: For each vr or 77 decay mode, the number of events A'^, the Gaussian width on the mass 
distribution of those data events, cr, the values of S, (3 (from MC), and the difference ((5)— /3, (see 
text). Uncertainties shown are statistical. 



Channel 


AT 


a 


{8) 


/3 


{5)-f3 






(MeV) 


(keV) 


(keV) 


(keV) 


77 (ttS) 


420 


3.51 


-285±195 


-51±21 


-234±196 


77 (7r°) 


4692 


1.94 


74± 46 


128± 8 


-54± 47 


77 


11140 


1.96 


419± 27 


126± 5 


293± 27 




1278 


2.83 


384±102 


233±18 


151±104 




3137 


1.12 


257± 24 


5± 4 


252± 24 




770 


0.91 


377± 44 


38± 6 


339± 44 



TABLE IL For each 77 channel, systematic uncertainties in (in keV) from the listed sources 
(see text); where applicable the degree of variation of the source level is given ("Var"). The sources 
marked with an asterisk (*) are assumed to be fully correlated across all modes; others are assumed 
to be uncorrelated. The final column combines the uncertainties across all modes with the weights 
given in the text. 



Source 


Var 


77 


37rO 7r^ 


^TT 7r° tt"* 


"tt 7 


All 


Fit Window 




14 


52 


12 


22 


9 




34 keV 


27 


32 


25 


32 


27 




11 keV 


9 


16 


9 


10 


9 


Bias 


A/3 


42 


78 


2 


13 


9 


Backgrounds 




19 


67 


17 


27 


12 


Pj^± scale 


0.01% 


1 


4 


5 


1 


3 


scale 


0.6% 


13 


26 


3 


7 


3 


MC Modeling* 




46 


46 


46 


46 


46 


Sum 




74 


132 


57 


68 


57 
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FIG. 1: Distributions of 6 for two evaluations of M^o using vr*^— 5-77 decay, with the data represented 
by the points with error bars and the Gaussian fit overlaid. The solid line portion of the fit indicates 
the window used for the fit and the dashed portions its extension. The solid line histogram 
represents MC simulation of signal and backgrounds with M^^=M^P^ and M^*^=547.78 MeV. 
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FIG. 2: Distributions of 6 for tj decay modes as shown. Symbols are as defined in Fig. 1. 
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